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ABSTRACT 

The He II A1640 emission line has been suggested as a direct probe of Pop- 
ulation III (Pop III) stars at high-redshift, since it can arise from highly ener- 
getic ionizing photons associated with hot, metal free stars. We use the HST 
WFC3/F130N IR narrowband filter to probe He II A1640 emission in galaxy 
IOK-1 at z = 6.96. The sensitivity of this measurement is >5x deeper than 
for previous measurements. From this deep narrowband imaging, combined with 
broadband observations in the F125W and F160W filters, we find the He II 
flux to be 1.2 ± 1.0 X 10"^*^ ergs s~^ cm~^, corresponding to a la upper limit 
on the Pop HI star formation rate (SFR) of ~ 0.5 Mq yr~^ for the case of a 
Salpeter IMF with 50 < M/Mq < 1000 and mass loss. Given that the broad- 
band measurements can be fit with a UV continuum spectral flux density of 
~ 4.85 X 10"^*^ X A~^'^^ ergs s~^ cm~^ A~^, which corresponds to an overall SFR 
of ~ 16l|6 ^0 y^~^) massive Pop HI stars represent < 6% of the total star 
formation. This measurement places the strongest limit yet on metal-free star 
formation at high redshift, although the exact conversion from He II luminosity 
to Pop HI SFR is highly uncertain due to the unknown IMF, stellar evolution, 
and photoionization effects. Although we have not detected He II A1640 at more 
than the 1.2a level, our work suggests that a > 3cr level detection is possible with 
JWST. 
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Introduction 



The first stars, i.e., Population III stars with zero metallicity, are likely to be very 



mass i ye, due to inefficient cooling and the large Jeans mass in the early Universe ( lAbel et al. 



2OO2I : lYoshida et al.ll2006l ). Massive Pop III stars are probably key contributors both of the 
ionization photons that caused c osmic reionization and to the early phase s of cosmic chemical 
evolution (IWhalen et al.l I2OO8I : iMackey et al.l 120031 : iBromm et al.l 120031 ). Although recent 
obser vations have provided major breakthr oughs in understanding the universe at 2; > 7 



(e.g., iBouwens et al.ll2010al : lYan et al .1120 111 ), detection of high redshift galaxies that contain 
Pop III stars remains a major challenge obser vationally. Indirect constrain ts on Pop III stars 
come from mapping the reionization history ( iFan. Carilli. &: Keatingll2006l ) or from studying 
the chemical evolution and the enrichment hi story through the observation of Extremel y 
Metal- Poor stars and Hyper Metal-Poor stars flBeers fc Christliebll2005t iBessel et al.ll2004l ). 
There are still no direct dectections of high redshift galaxies hosting star formation in zero 
metallicity gas. 

Models show that Pop III stars are considerably hot ter than present-day stars, for a 
given stellar mass (e.g., iJohnsoru I2OIOI : iBromm et al.l 120031 ). Emission of ionizing radiation, 
specifically photons with energies above 54.4eV that ionize He II, is enhanced by the high 
surface temperature photospher es of Pop HI stars. He II A1640 ha s been long sug gested as 
a direct probe of Pop III stars jTumlinson fc Shulll I2OO0I : lohlboOli ISchaereilEooi ). Hell A 
1640 emission is also relatively easy to model compared to resonance lines such as Lya or He 
II A304. The optical depth to He II A1640 is small in almost all circumstances, so it does not 
suffer the radiati ve transfer effe cts that complicate the interpretation of Lya flux in high- 
redshift galaxies flSchaererll2003l). In addition, it does not suffer from IGM Gunn- Peterson 
absorption as Lya does ( jTumlinson fc Shulll l2000l ). 



How can He II A1640 emission from star formation in zero metallicity gas be distin- 
guished from other sources ? Calculations show that the He II A 1640 emission in a metal- 
poor stellar population with metallicity as low a s Z~10~^ is still <10~ ^ times lower than that 



in a population formed in zero metallicity gas ( jTumlinson et al.ll2003t ISchaererl |2003| ). It is 
also known t hat nebular He II A4686 was detected in a fairly large fraction of metal-poor HII 
regions (e.g. , ISkillman fc Kennicuttlll993[ ) . The follow-up observations and theoretical inves- 



tigations suggest that Wolf-Rayet (W-R) stars are responsible for this nebular He II emission 
( Ide Mello et al.lll998l ). but the hardness of this emission more than one order of magnitude 
weaker than that originated by massive Pop HI stars (ISchaererl l2003t I Jimenez fc Haiman 
20061 ). Also, the He II A1640 equivale nt width of massive Pop HI stars could be a factor of 
2 larger than that of the typical AGN dPrescott et aubood : IeIvIs et al.lll994h . 



A number of studies have been carried out to search for He II emission in high-redshift 
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galaxies at 4 < z < 6.6. IPawson et al.l (120041 ) . lOuchi et all fcOOSf ). and iNagao et all fl2005[ ) 
searched for He II A1640 in either stacked or individual spectra of Lya emitting galaxies. 
Nagao et al.l (|2008[ ) carried out a survey for Lya-He II emitters using a combination of inter- 
mediate and narrow-band filters in the optical window. But these observations so far have 
yielded only nondetections, constraining the massive Pop III star formation rate (hereafter 
SFRpnpTTi) to a few Mf7) yr ~^, usually a few tenths of the overall SFR in these galaxies. Yet 
Jimenez fc HaimanI ( l2006l ) predict SFRpopm still to be signficant at z ~ 3-5. Even stronger 
metal-free star formation is expected at higher redshifts, due to less chemical feedback from 
Pop III stars after their initial burst. However, at z > 7, He II is in the J-band, mak- 
ing ground-based observations extremely challenging. Due to its large field of view, high 
throughput, and low background, the new IR channel of the Wide Field Camera 3 (WFC 
3) on HST is the most powerful tool for detecting galaxies at z> 7 and for studying galaxy 
evolution at the end of cosmic reionization. 

In this paper, we report the strongest upper limit yet on the contribution of Pop HI stars 
in high-redshift galaxies, using deep narrow band imaging of galaxy IOK-1 ( llye et al.l 120061 ) 
at z = 6.96 in the HST/WFC3 F130N band, which includes He II A1640 at this redshift. We 
have also carried out deep observations in broad bands to measure the continuum level and 
UV-based total SFR. In §2, we discuss the observations and data reduction. Photometric 
results, including the limit on He II flux, are presented in §3. We discuss the implication 
of our results on Pop HI star formation in IOK-1 in §4. Throughout this paper, we adopt 
a cosmology ba sed on the fifth year Wilkinson Microwave Anisotropy Probe (WMAP) data 
komatsJboOoi ): = 0.72, Q„^ = 0.28, Qb = 0.046, and Hq = 70 km s"^ Mpc'^ 



2. Observation and the data reduction 

Our target, IOK-1, was discovered in a narrow band Lya emitter survey using the NB973 
filter on the Subaru Suprime-Cam (lye et al. 2006), which is centered on Lya at 2; ^ 7.02. 
This galaxy is spectroscopically comfirmed at z = 6.96. IOK-1 has a total Lya luminosity 
of 1.1 ± 0.2 X 10^2 erg s-^ corresponding to a Lya-based SFR SFRLya ~ 10 ± 2 Mq yr-\ 
(lye et al. 2006). Q 

In the observed frame, the He II A1640 emitted by IOK-1 is located at 13054 A, which 
is in the passband of the WFC3 F130N narrow band filter (Aq = 13006 A and FWHM = 



^These error bars only reflect the photonietrie uncertainties. There are big systematic uncertainties due 
to Lya radiative transfer effect, galaxy IMF and metallicity assumptions. The SFRLya could change by a 
factor of 2 given these systematic uncertainties (also see §4). 



-4- 



150 A). Resonant scattering of Lya photons and IGM absorption of the blue wing usually 
systei natically redshift th e Lya line center up to ~ 700 km s~^ against the systemic velocity 



(e.g., 



Steidel et al.ll2010l ). Given that He II will be located slightly blueward of the peak 
inferred from Lya peak and have a narrower line width, we expect that the entire He II 
line falls securely into the most sensitive part of the F130N filter. IOK-1 was observed by 
HST WFC3 in March 2010. Eight orbits (~ 20, 000 sec integration) were devoted to the 
F130N filter to measure the He II flux; we also observed this galaxy using the F125W filter 
for two orbits to determine the underlying continuum level and UV continuum-based SFR. 
A F160W band observation was also carried out (Egami et al. 2011, in preparation). 

The F130N observations are divided into two visits, each consisting of a 4-point dither 
sequence with one-orbit per dither. During the the first visit, the first two dither positions 
were affected by the presence of a bright ghost image very close to IOK-1. We exclude these 
two images from further analysis. The F125W continuum observation was a two-orbit single 
visit, which is divided into two identical 4-point dither sequences, one per orbit. No ghost 
image affected the F125W observations. The individual images were reduced by WFC3-IR 
standard pipelines. Bot h the F130N and F125W images were combined using Multidrizzle 



(IKoekemoer et al.ll2002l ) with final^cale = 0'.'06, which is 0.48 of the original pixel size, and 
finaLpixfrac = 0.7. High resolution images from the F130N and F125W bands are shown in 
Figure 1. 



3. Results 
3.1. Photometry 



Photometry is performed with SExtractor (IBertin fc Arnoutslll996l) using the rrns ma p 
converted from the inverse variance image generated by Multidrizzle (jCasertano et al.ll2000l ). 
We measure the fluxes in the broadband (F125W and F160W) and the narrowband (F130N) 
images using the same Kron-like elliptical aperture determined from the F130N image (black 
elliptical aperture in Figure 1). The results are listed in Tabled] We also measure the flux 
in the F160W image of IOK-1 (Egami et al. 2011, in preparation), using the same aperture 
as used in the F125W and F130N images, flnding a flux density of /pieow^ 2.44 ±0.14 x 
]^g-2o gj,g g-i corresponding to a magnitude of mpieow = 25.98 ±0.06. We fit the 

photometry by assuming a model spectrum with a power law continuum with a dimensionless 
flux density fcon = ^(-^/lA)'', where A is a constant, and a narrow Gaussian for the He II 
emission line of flux Fneii. The best-flt SED and photometry are shown in Figure 2. We 
flnd: 



/con(A) = (4.85 ± 0.23) X 10-^°(A/1 A)" 



2.46±0.36 
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Fhcii = (1.2 ±1.0) X 10 



_ 1 Q 1 9 

ergs s cm . 



(2) 



k.t z = 6.96, this corresponds to a total Hell luminosity of ^Heii = Fhcii x 47rDL = (6.6 
±5.5) X lO^^ergs s~^, or a 2-cr upper limit of 1.11 x 10^^ erg s~^. This 2-cr upper limit is a factor 
of >4-5 X deeper than previous measurements. The rest-frame equi valent width of He II is 
4.2±3.5A.| Our results also show that IOK-1 has a blue continuum. iBouwens et al.l (j2010b( ) 
studied the value of the UV- continuum slope (3 in the Hubble Ultra Deep Field (HUDF). 
For luminous L*,=^ galaxies f lSteidel et al.lll999[ ). /3 ~ —2.0 ± 0.2. For lower luminosity 

0.1L*^3 galaxies, /3 3.0 ± 0.2. IOK-1 has an absolute magnitude Mab(1500A) = -21.3, 

comparable to L*,^^. 



3.2. Morphology 



The half-light radius ri/2 of IOK-1, based on our SExtractor measurements flBertin fc Arnouts 



19961 ) and corrected for WFC3 PSF broadening, is 0'.'12 in F125W. This corresponds to 

0. 62 ± 0.04 kpc. IOK-1 is an extremly compa ct galaxy, consiste nt with the size of L*^^^ 
galaxies at z ~ 7 detected in the HUDF09 (jOesch et al.l |2010| ). The observed surface 
brightness of IOK-1 is /ij ~ 24.5 mag arcsec"^, corresponding to a surface brightness of 
^rest ~ 15.5 mag arcsec"^, after (l±z)^ correction for cosmological dimmi ng. From Figure 

1, it is obvious that IOK-1 consists of two components. We use GALFIT ( iPeng et al.l 120021 ) 
to deblend the F125W continuum image, assuming an exponential profile for both compo- 
nents. The F125W image shows roughly equal brightness for the two components. The north 
west component has an effective radius of 0.49 ±0.04 kpc; and the southeast component has 
an effective radius of 0.57 ± 0.03 kpc. The two components are projected about 0'.'2 (~1 
kpc) away from each other. More detailed discussions on the morphology are presented in 



Note that our observations could only detect He II flux stronger than 3xl0~ erg s 
level, corresponding to a He II A1640 equivalent width of '^llA. which is a factor of 6 larger than tha t 
of normal stellar populations in the Lyman break spectrum of stacked z^3 galaxies (jShaplev et al.l 120031 ). 
and a factor of > 3 grea ter than the simulated maximum He II equivalent width generated by W-R stars 
(jBrinchmann et al.ll2G08[ ). 



Table 1: SExtractor photometry in different bands 



Filter 


Flux Density f^ 


niAB 




(10"2° erg s^^ cm~2 A"^) 




F130N 


4.39 ±0.68 


25.42 ±0.17 


F125W 


4.07 ±0.19 


25.59 ±0.05 
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Egami et al. (2011, in preparation). We also deblend the F130N narrow band image, which 
contains He II A1640 emission. The north-west component is somewhat brighter. How- 
ever, given the low signal-to-noise ratio of the data, this difference, while intriguing, is not 
statistically significant. 



Discussion 



What limit can our observations place on the SFRpopin and its contribution to the 
overall star formation in IOK-1 at z ~ 7 ? Photoionizatio n by AGN or hot dense s tellar 



winds from Wolf-Rayet stars (see e.g., IShapley et al.l l2003t iBrinchmann et al.l |2008[ ) may 



contribute to He II A 1640 emission although at a lower level than metal-free stars. More 
importantly, massive stars with the lowest, but non-zero metallici ty (< 10~^) cou ld generate 
a comparable amount of He II A1640 emission to Pop HI stars (jSchaererll2003l ). However, 
here we assume that the observed He II emission originates entirely from metal-free stars. 
Therefore, our derived SFRpopm should be regarded as an upper limit. In the case of constant 
star forma tion, at equilib rium, recombination line luminosities Li are proportional to the 
SFRpopiii JSchaereiibooi ). so 



^Hell 



Cl640(l - /esc)Q(He^ 



( SFRpopiii \ 
V Moyr-i ) 



^1640, norm 



/ SFRpopIII \ 

V Moyr-1 ) 



(3) 



where Ivi64o,nomi = Ci64o(l — /esc)Q(He"'~) is the theoretical He II A1640 line luminosity 
norr nalised to SFR= 1 yr~^. Ci64o is the He II A1640 emission coefficient given in Table 
1 of ISchaererl (120031 ): ci64o = 5.67 x 10~^^ ergs for Te = 30fcK, where nebular emission is 
calculated assuming case B recombination. (^(He"*") is the number of He^ ionizing photons 
per second, ff.sc represents the fraction of total ionizing radiation released into the IGM 
without being coupled to the ISM in the g alaxy. We assume f ^. y = 0, as it is e xpected to 



be small among high-redshift galaxies (e.g.. iGnedin et al.ll2008l ). ISchaererl (120021 ) calculates 
the total He"*" ionizing photon flux for a Salpeter IMF with a range of lower and upper mass 
cut-offs and different assumptions of mass loss. 

The conversion from He II A1640 luminosity (Lhoii) to SFRpopm is affected by a num- 
ber of uncertainties, including: (1) IMF, (2) mass loss, and (3) photo-ionization model. 
For the IMF, theoretical calculations show that the lack of efficient cooling mechanisms 



generally result in an extrernely top-heavy IMF compared to a galactic IMF (lAbel et al. 



2002; Bromm fc Larson 2004; Yoshida et al. 2006; O'Shea fc Norman 2006). However, both 



the shape and mass range of such an IMF are poorly constrained. IScannapieco et al.l (120061 ) 
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Table 2: Pop III star formation rate 



IMF(Salpeter) 


Mass Loss 


-^1640,nomi 

(ergs S^^) 


SFRpopiii 

(Mo yr-i) 


1 < M/Mq < 500 


No 


9.66 X 10^° 


6.8 ±4.7 


50 < M/Mq < 500 


No 


6.01 X 10^1 


1.1 ±0.9 


1 < M/Mo < 500 


Yes 


3.12 X 10^1 


2.1 ± 1.5 


50 < M/Mq < 1000 


Yes 


2.33 X 10^2 


0.3 ±0.2 



suggests a lower mass cutoff of 0.8 based on metal abundances in Galactic halo stars. Nu- 
merical simulations su ggest that stars as mass ive as 500 M0 can be formed by accretion into 
a primodial protostar fjBromm fc Larsonll2004r), which is gene r ally taken as an upper limit in 



differ ent simulations ( ISchaerei 



2OO2I: 



Scannapieco et al 



20031 



Tum 



2010[) . Both a Salpet er IMF dSchaereil 12003 ! IScannapieco et aP I2OO3I ) and a Log-normal 



linson 2006; iRaiter et al. 



IMF ( iTumlinsonl l2006l ) have been used in model calculations. Table 2 shows a factor of six 
difference in inferred SFR arising from different lower mass cutoffs in the Salpeter IMF. 

As shown in Table 2, models with strong mass loss have a larger conversi on factor 
L^f\4( j r^nr•m', therefore, for th e same Hell luminosity, they yield a lower SFRpopm ( ISchaerer 



2OO3I ). iKudritzkil (j2002[ ) shows that very low-metallic it y sta rs close to the Eddington 



limit are subject to non-zero mass loss. ISmith fc Owockil ( l2006l ) argue that massive shells 
around luminous blue variables and so-called supernova impostors indicate that continuum- 
driven winds or hydrodynamic explosions dominate the mass loss of very massive stars, 
a claim which is insensitive to metallicity and so could apply to Pop 111 stars. Ekstrom 
(2007) explores the effects of rotation, anisotropic mass loss, and magnetic fields on the 
core size of a Pop 111 star, pointing out that, under certain conditions, very massive stars 
(I4OM0 < M < 26OM0) losing mass through rotation could avoid ending their lives as 
pair- instability supernovae (PISN). Whether Pop 111 stars end their lives as PISN or not will 
affect the metal ejection efficiency and hence affect the SFRpopm. From Table 121 stronger 
mass loss will lower SFRpopm by a factor of more than three. 



The derived SFR also depends on the photoionization model. iSchaererl (120021 ) cal- 
culates nebular emission lines using standard case B recombination, in which the optical 
depth of the Hll region is large. Eq. (3) is derived under this assumption. However, for 
low metallicity nebulae ionized by very hot Pop 111 stars, the case B predictions for line 
and continuum era ission may have non-negligible deviations from real nebular astrophysics 
( IRaiter et al.ll2010l ). Therefore, the SFRpopm upper limits d erived from Eq. (3) then need 
to be qualitatively revised here. According to IRaiter et al.l ( 120101 ). for a given He 11 lumi- 



nosity, SFRpopiii could be a few times higher, depending on the ionizing parameters, e.g.. 
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the hydrogen number density n{H), inner radius of the nebula 
photon flux Q(H). 



and hydrogen ionizing 



Our upper limit on the He II emission coombined with our detection of IOK-1 in the 
Spitzer IRAC 1 and 2 bands (Egami et al., in prep.) suggest that IOK-1 is dominated by 
stars of metallicity above the critical value (Zcrit ~ 5 x lO'^^ZpX lead ing to a normal IMF 



(e.g., Bromm 2004). Using the relation given in iMadau et al.l (119981 ) . the rest frame UV 
luminosity of the galaxy IOK-1 corresponds to an overall SFR ~ 20to 9 y^~^- These 
error bars only reflect the uncertainties in the J-band magnitude. The true errors are larger 
as the conversion from UV luminosity to overall SFR is highly and systematically uncertain. 
In this conversion, Madau (1998) assur aes solar metal licity; lower metallicities down to 0.0004 
could lower the SFR by a factor of 2 (jSchaererll2003l ). Furthermore, the conversion assumes 
a Salpeter IMF with mass ranging from 0.1 to 125M0, whereas assuming a Scalo IMF would 
double the derived overall SFR. The UV-based overall SFR which should be less obscured 
is som ewhat higher than that based on the Lya emission line (~ 10 ± 2 M© yr"^; live et al.l 
20061 ). but consistent given the uncertainties in both the stellar population models and the 
Lya radiative transfer effects. 

The conversion of UV luminosity to an overall SFR above assumes a conversion factor 
between UV luminosity and SFR for Pop I and Pop II stellar populations, without correcting 
for Pop III stars with a top-heavy IMF. Using the fiducial Pop III Salpeter IMF of 50 - 500 
Mq and no mass loss, we find that our derived upper limit of 2.0 Mq yr~^ for SFRpopm 
corresponds to an upper limit of 1.2 x 10"^" ergs s~^cm' -2 A"^ in F125W flux density. Note 
that the observed flux density is 4.07 ±0.19 x 10"^'' ergs s~^cm~^ A~^. After correcting for 
the UV-luminosity contributed by Pop III stars, the overall SFR is ~ 16l|g Mq yr~^ ( As 
stated previously, this result could change a factor of 2 given the systematic uncertainties 
associated with the IMF and metallicity assumptions). Thus, IOK-1 is not dominated by 
very massive Pop III stars. It should be noted that low mass metal-free stars do not produce 
an appreciable amount of He II A1640 emission. Therefore, our observations only constrain 
the amount of high-mass Pop HI star formation in the galaxy and cannot be used to rule 
out the existence of low mass, low surface temperature Pop HI stars (alth ough, as discussed 
in_§4, most theoretical calculations favor a top heavy IMF for Pop HI). iTrenti fc Stiavelli 
(120091 ) show that the fraction of Pop HI s tars /ttt(^) per da rk matter halo at 2; ~ 10 is only a 
few thousandths. iDave et al.l (120061 ) and iTumlinsonI (120061 ) also suggest that fiii{z) in halos 
at 2; ~ 7 is small: < 1%. Taking into account the systemic uncertainties in the UV-based 
overall SFR, our observations indicate that the galaxy IOK-1 cannot be dominated by Pop 
HI stars with very top heavy IMFs. For example, for a Salpeter IMF with 50 < M/Mq < 500 
and no mass loss, the ratio of SFRpopin to the overall SFR is < 25%. For a Salpeter IMF 
with 50 < M/Mq < 1000 and mass loss, the ratio is < 6%. 
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Our deep HST narrow-band imaging places the strongest constraint yet on the SFRpopm 
in high-redshift galaxies. Although we have not detected the He II emission line in IOK-1 
at more than 1.2a, this limit suggests that detection of He II emission from Pop HI stars 
may be already within the reach of current facilities. Future facilities, especially JWST and 
GSMT, will be likely to probe down to the > 3a level and directly detect the signatures of 
the earliest star formation in the universe. 
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-VI- 




Fig. 1. — High resolution images of galaxy IOK-1 in the F125W (left) and the F130N (right) 
bands, with contours spaced by 1.4 sky rms for the F125W image and 1 sky rms for the 
F130N image. The black elliptical aperture is determined by the F130N image. IOK-1 is 
clearly resolved into two components separated by ~ 0'.'2. Photometric analysis of the F130N 
image reveals a flux density excess of 1.2 a. 
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Fig. 2. — Best-fit spectrum (black line) of galaxy IOK-1 with the total UV continuum, as well 
as the Lja and He II emission lines. The filter response curves of F125W (blue dashed line), 
F160W (dark blue dashed line) and F130N (brown dashed line) are plotted. In addition, 
photometry in three different bands is overplotted at the effective wavelength of each filter. 



